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Photoabsorption Spectra of Na+ clusters: Thermal Line-Broadening Mechanisms 
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Photoabsorption cross sections of small sodium cluster cations (Na+, n=3,5,7 and 9) were cal- 
culated at various temperatures with the time-dependent local-density-approximation (TDLDA) 
in conjunction with ab initio molecular dynamics simulations, yielding spectra that agree with 
measured ones without ad-hoc line broadening or renormalization. Three thermal line-broadening 
mechanisms are revealed: (I) lifting of level degeneracies caused by symmetry-breaking ionic mo- 
tions, (II) oscillatory shifts of the entire spectrum caused by breathing vibrations, and (III) cluster 
structural isomerizations. 
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Optical spectroscopy provides invaluable insights into 
the electronic structure, ionic configurations, thermal 
processes and dynamics in metal clusters, as well as 
about the size-dependent evolution of these properties. 
Indeed, investigations of these issues have been pursued 
rather intensively for over a decad&,both experimental! 
with photodepletion spectroscopyicl and theoretica 
with the main methodology employing thc-TJDLD 
in conjunction with cither jcllium modclsutl (including 
allowance for shape-deformations) where the ionic back- 
ground is smeared out uniformly, or with electronic struc- 
ture calculations of various degree of sophistication where 
the discrete nature of the ions is incorporated accu- 
ratelyLm or perturbativelyB. While valuable information 
pertaining to the optical excitations and damping mech- 
anisms in metal clusters has been obtained through such 
studies, including the sensitivity of optical features to 
cluster geometries, a first-principles theoretical descrip- 
tion of the optical line-shapes, the absolute absorption 
cross-sections, the relevant microscopic line-broadening 
mechanisms and their thermal variations, is lacking. This 
is reflected in the common ad-hoc employment of line 
broadening through convolution of spectral lines calcu- 
lated for selected static clusteLcanfigurations with Gaus- 
sian or Lorentzian functionsETQ'Lru, or (multiplicative) 
renormalization of the calculated spectra to the measured 
onesBEl. 

In this letter, we demonstrate that photoabsorption 
cross sections, (t(oS), calculated via the TDLDA along 
finite-temperature Born-Oppcnhcimcr (BO) local-spin- 
density (LSD) molecular dynamics (MD) (BO-LSD- 
MDEl) phase -space trajectories, provide a quantitative 
ab-initio description of the absolute magnitudes, peak po- 
sitions, and line shapes of optical absorption spectra mea- 
sured at various temperatures from Na+ (n=3, 5, 7 and 
9) clustersa. Furthermore, analyses of the results, and 
in particular the temporal correlations between the cal- 
culated ^(lu) (where the superscript denotes the value 
at instant t) and the dynamically evolving cluster con- 
figurations, reveal three thermally induced mechanisms 



which govern the spectral characteristics: (I) symmetry 
breaking of the optimal (0 K) cluster geometries due to 
thermal motion of the ions, leading to lifting of electronic 
level degeneracies and resulting in splitting (fragmenta- 
tion) of spectral lines; (II) symmetry conserving breath- 
ing motions that modulate the effective density, result- 
ing in oscillatory frequency shifts of the entire spectrum 
("spectral sweeping") that influence the optical absorp- 
tion peak positions and line width; (III) opening of the 
accessible configurational space, resulting in structural 
isomerization and consequent electronic structure and 
spectral variations. The first two mechanisms are opera- 
tive for smaller clusters (Na^~ and Na^) even at relatively 
low temperatures (e.g. 100 K), while all three occur for 
the larger ones. Furthermore, the thermal dependence 
of the spectral characteristics is analyzed, resulting in 
remarkable first-principles agreement between the calcu- 
lated and measured optical absorption cross sections. 

The Kohn-Sham (KS) .equations with generalized gra- 
dienJ—corrections (GGA)li3 and non-local pseudopoten- 
tialaiJ were solved for Na+ clusters (n=3, 5, 7, and 9), 
and based on the Hellman-Feynman forces, finite temper- 
ature trajectories were generated by Langevin dynam- 
ics through the use of the BO-LSD-MD (with GGA) 
method! 2 ]. For each cluster size, starting from the (pre- 
sumed!) GGA-optimized ground state (GS) geometry 
and an equilibration of 10 ps duration, the system was 
allowed to evolve dynamically at the prescribed tempera- 
ture for an additional period of i s =10 ps. At each instant 
t, the excitation energies lo\ were, determined by solving 
the TDLDA eigenvalue problerr£3 
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where i, j and k, I denote particle-hole pairsO, and e*^ the 
corresponding KS orbital energy differences. The cou- 
pling matrix is given by 
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with the GS electron density p l , the product of two (real) 
KS orbitals <p\j , and the LDA exchange-correlation po- 
tential v xc . The oscillator strength (OS) 
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corresponding to the transition energies £j|, were accu- 
mulated in bins of width Aw=0.025 eV yielding the ab- 
solute photoabsorption cross section per valence electron 
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We note first the overall good agreement between the 
magnitudes and shapes of the calculated and measured 
photoabsorption cross sections (shown here for Na+, 
n=3,5 and 9); compare the red histograms with the black 
curves in Fig. 1, and note that neither (artificial) numer- 
ical broadening nor a renormalization of the spectra has 
been applied. The slight deviations in peak positions ace 
within the known accuracy of the TDLDA (±0.2 eV)EH. 

The calculated Na^j" spectrum at TsslOO K (Fig. la) 
exhibits a low energy peak originating^ from two transi- 
tions a4=w|«2.65 eV, from the occupied s-like orbital to 
two empty p-like KS orbitals located in the plane of the 
Nag equilateral triangle (see grey insets 1, 2 and 3 at the 
left of Fig. la). The high energy peak at ^3=3. 41 eV is 
due to an excitation to another p-like orbital that is per- 
pendicular to the cluster plane (grey inset 4 in Fig. la). 
In the experiments cr(uj) is determined from the deple- 
tion of the Na+ intensity due to dissociation following 
absorption of a photon hto. Excitation into the in-plane 
anti-bonding orbitals (insets 2 and 3) promotes such dis- 
sociation and is accompanied by exhaustion of the OS for 
the first (lower energy) measured peak. Excitation into 
the out-of-plane orbital (inset 4) has no direct destabi- 
lizing effect and consequently in the measurements .only 
66% of the Thomas-Reiche-Kuhn sum rule (TRK)ll9 is 
observed while the full TRK sum rule is found in the 
theoretical spectrum shown in Fig. la (as well as in all 
the calculated spectra that we show here). 

Thermal motions distort the T>3h symmetry of the Najj" 
GS and the degeneracy of the low-energy transitions lo\ 
and ui2 is lifted. Indeed, the evolution of 5u>2i=oj\—oj\ 
coincides with the temporal behavior of the standard de- 
viation Ak=^/^j(&* — 6*) 2 /3 of the three bond distances 

b\ (see Fig. 2a). Additionally, the mean bond distance 
b is strongly anti-correlated with the mean frequency 
c^i2 = ( w i +ujI)/2 as well as with the position of u)\ (com- 
pare the solid line with the lower and upper dashed curves 
in Fig. 2b) - that is, an increase of the ionic density (de- 
crease of &') results in a blue shift of the entire spectrum. 
Thus, the spectral line width of Nag" can be fully un- 
derstood in terms of line-splitting and spectral sweeping 
with the first caused by symmetry breaking (degeneracy 
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1. Theoretical (red histograms) and experimentalEI 



(black curves) photo absorption cross section for NaJ at 100 
K (a), NaJ at 100 K (b) and 300 K (c), as well as for Na<[ 
at 100 K (d) and at 450 K (e). Grey insets in (a) display 
isosurfaces of the 4 lowest KS orbitals where yellow and pur- 
ple distinguish the sign of the wave function and light blue 
spheres mark the position of the ions. The right top cor- 
ner insets in (b-e) show static TDLDA spectra: (b) for the 
NaJ GS (structure a), and for an instantaneous structure (3 
with an elongated left triangle; (c) for an instantaneous bent 
planar geometry (7). The calculated spectrum shown in (d) 
derives exclusively from configurations lying in the basin of 
the GS structure a whose oblate spheroidal electron density 
is depicted by the purple isosurface. The static spectrum of 
a is shown at the top inset in (d), and in the bottom inset 
we show the decomposition of that spectrum into the com- 
ponent along the z-axis (green sticks) and the perpendicu- 
lar components (red sticks). The Nag" isomer (3 is shown in 
(e), along with its electron density (purple isosurface) and its 
static spectrum (right top inset). To highlight the sensitivity 
of the experimental results rto the temperature, we include in 
(e) the measured spectrumB at 540 K (grey curve). 
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lifting, mechanism I) and the latter by breathing vibra- 
tions of the cluster (density oscillations, mechanism II). 

Similarly, the finite temperature dynamics of Na^~ dis- 
torts (even at T=100 K) the optimal D2d GS symme- 
try (see GS a in Fig. lb where the perpendicular left 
and right triangles are equivalent, and compare to the 
instantaneous structure (3 shown in Fig. lb with an elon- 
gated left triangle) , resulting in fragmentation (I) of the 
GS line at uj= 2.8 eV into two separate spectral lines ui u 
and ujd Ell (compare red and blue lines in the up-right 
inset of Fig. 2b). Consequently, the temporal variation 
of the absolute difference of the average bond lengths in 
the left and the right triangle (A* r = 



\b\- 



serving as a 



measure of the thermal structural distortion) correlates 
directly with the spacing between the fre- 

quencies of the spectral fragments (see solid and dashed 
curve in Fig. 2c). 

While the relatively moderate structural distortions 
occurring at 100 K are not sufficient to lift the degeneracy 
of the states associated with the low energy absorption 
line of Na^" (at ~2 eV, compare insets a and (3 at the 
upper right of Fig. lb) which remains relatively sharp, 
increase of the temperature to 300 K results in bent con- 
figurations (see instantaneous structure 7 on the left in 
Fig. lc) where the low energy line is fragmented (see 
inset 7 at the upper right of Fig. lc), resulting in en- 
hanced broadening of that spectral line (compare Figs, 
lb and lc). Further broadening of the lines is caused by 
the aforementioned breathing modes (spectral sweeping 
mechanism II) reflected in the direct correlation between 
the interionic Coulomb repulsion E l coul (as a measure 
of the ionic density) and the average excitation energy 

^*=Ej /M/(« - !) shown in Fi g 2d for the T = 300 K 
simulation. 

For larger clusters, the occurrence of thermal iso- 
merizations opens an additional line-broadeiiing channel 
(III)El The GS structure of the Na^ clusterEJ is shown at 
the left of Fig. Id (marked a) and the TDLDA spectrum 
calculated for this static configuration exhibits several 
absorption lines distributed in a bimodal-like manner (see 
top inset at the upper right of Fig. Id). This bimodality, 
which is in contradiction to-the single-line spectrum pre- 
dicted by the jellium modeH, originates from the "oblate 
shape" of the cluster with radii i? 2 =2.68 A along the D$h 
symmetry axis (z in Fig. Id) and R xy =3.22 A in the xy 
plane (these radii were determined from the diagonalized 
moment of inertia of the ionic GS structure a). The to- 
tal electron density (see purple isosurface in Fig. Id) is 
almost an ideal oblate spheroid. A decomposition of the 
OS of the T>3h GS into a z (^=3 term in eq. |l|) and the xy 
contribution (v=l and 2 terms) shows clearly that, as ex- 
pectedu, the smaller spatial dimension in the z direction 
is associated with a higher excitation energy (compare 
green and red line in bottom inset at the upper right of 
Fig. Id). The bimodal character of the 100 K photoab- 
sorption spectrum is evident in both the measured and 
calculated spectra (Fig. Id), with further broadening 
caused by thermal motions through line-fragmentation 
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FIG. 2. Evolution of structural and spectral quantities with 
time: (a) the standard deviation of the three Najj" bonds 
(solid curve) coinciding with 60)21 (dashed curve) at 100 K; 
(b) the mean bond length 6* of NaJ (solid curve) and the 
mean frequency a)* 2 (lower dashed curve), as well as the tran- 
sition energy uj 3 (upper dashed curve) at 100 K; (c) absolute 
difference A, r =|6 ; — 6*| between the average bond length in 
the left and the right triangle of Na^ (solid curve) and Sui^ 
(dashed curve) at 100 K; (d) Interionic Coulomb repulsion per 
ion in Na^~ (solid line) and mean energy of the corresponding 
TDLDA spectrum w* (dashed curve) at 300 K; (e) A^ 5 (solid 
line) and r* (dashed line) for Na+ at 450 K. 



and breathing vibrations of the cluster (mechanisms I 
and II). 

The main thermal effect observed in the spectra of Na^J" 
is the conversion of the low temperature (e.g. T=100 K, 
Fig. Id) bimodal spectrum to one with a single broad 
maximum at higher temperatures (e.g. 450 K, Fig. le). 
This change in the spectrum is caused mainly by trans- 
formations between the GS geometry (a, Fig. Id) and 
the structural isomer [3 (see configuration on the left 
of Fig. le). The static spectrum of the (3 isomer is 
shifted to lower energies with respect to that of the GS 
(a) structure (compare top insets marked a and (3 at 
the upper right corners of Figs. Id and le). Thus, the 
high-temperature broad spectral feature (Fig. le) is due 
to contributions from both the a and (3 isomers to the 
phase-space trajectories, generated by the Langevin MD 
simulation at 450 KBEL 

The /?-isomer may be characterized as having an oc- 
tupolar shape, reflected in the shape of the electron den- 
sity shown in Fig. le (see purple isosurface). To quanti- 
tatively analyze the spectral broadening caused by tran- 
sitions (at high T) between the GS (a) and the /3-isomer 
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(which involve changes in shape-multipolarity) we cal- 
culate the dimensionless multipole shape parametersE3, 

ahn= 3r t^f 3 +i Y%=i r i Y tm(Oi,<t>i) where (n,Qi,<f>i) are 
the polar coordinates of the ith ion with respect to 
the center-of-mass of the cluster and r s is the den- 
sity parameter of bulk Na. We defined a parameter 

Ae =Y^z'=2 T^m= -v l°^'-m| 2 that measures the deviation 
from sphericity due to the multipole modes up to the or- 
der I. Eig. 2e shows a remarkable correlation (i? 2 =84%, 
see ref.E3) between A\ =b and the width T* of the spec- 
trum at time t (T f is twice the standard deviation of 
a (w) for 2 < lu < 3.5 eV). Furthermore, including in 
such analysis only quadrupole and octupole deformations 
(i.e. A\ =3 ) results in a similar high value of i? 2 (73 %), 
while limiting consideration to the quadrupole deforma- 
tion alone (i.e. A\ =2 ) yields a poor correlation (i? 2 =19 
%) pointing to the strong influence of octupolar defor- 
mations on the optical response of metal clusters (even 
closed shell ones i.e. Nag"). This analysis confirms our 
conclusion pertaining to the importance of the isomeriza- 
tion mechanism (III) in explaining the thermal evolution 
of photoabsorption spectra (particularly at high temper- 
atures) . 

In summary, we demonstrated that calculations of ab- 
solute photoabsorption cross sections using the TDLDA 
in conjunction with ab initio MD simulations allow first- 
principles quantitative description and interpretation of 
optical spectra measured for sodium cations at various 
temperatures. Spectral line fragmentation resulting from 
symmetry-breaking ionic motions that lift electronic level 
degeneracies, frequency shifts (spectral sweeps) of the en- 
tire spectrum due to symmetry-conserving breathing vi- 
brations, and structural isomerizations are identified as 
the main thermal line-broadening mechanisms with the 
first two operative for the smaller cluster (Najj~ and Na^") 
and all three for larger ones. 
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